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complete multiplets. We derive the soft supersymmetry breaking terms and show that the gaugino 
sector can be parameterized by only two parameters. These parameters can be chosen appropriately 
to obtain various patterns of gaugino masses and different ratios among them. The sfermion sector 
can also be characterized by two independent parameters which can be adjusted to change the 
relative masses of squarks and sleptons. A judicious choice of parameters also allows us to achieve the 
lightest Higgs boson mass about 125 GeV. In this paper we focus on a scenario where a comparatively 
large hierarchy exists between the U(l) and SU(2) gaugino mass parameters. In such a case, the 
lightest Higgs boson originating from the decay of the next-to-lightest neutralino, following the 
direct production of chargino neutralino pair, can be considerably boosted. We show that a boosted 
Supersymmetric Higgs signal with a decent signal to background ratio can be obtained using the 
jet substructure technique at LHC with 8 TeV (14 TeV) center of mass energy and an integrated 
luminosity of about 30 fb _1 (15 fb _1 ). 
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I. INTRODUCTION 



The performance of the Large Hadron Collider (LHC) at CERN has been extraordinary in the last year and data 
of about 5 fb _1 has been already collected by each of the two major experimental groups CMS and ATLAS. With the 
data already accumulated, the LHC has left the LEP and Tevatron results far behind in many of the search channels. 
But, apart from the recent hint [1, 2] about a possible 124 - 126 GeV Higgs boson, no other signal of any new physics 
(NP) has been reported yet. 

The Minimally Supcrsymmetric Standard Model (MSSM) which is undoubtedly one of the most studied models 
beyond the Standard Model (SM) has been a benchmark for the LHC searches and limits on Supcrsymmetric particles 
have been shown by both the ATLAS and CMS collaborations [3] . 

Supersymmetry (SUSY) [4, 5] which is a beautiful theoretical idea by its own right has been in the spotlight 
among the phenomenologists for a long time because of its ability to tame the quadratic divergence in the scalar 
sector of the SM elegantly. However if SUSY has to be realized in nature, it must be broken so that it does not 
contradict experimental observations. From a theoretical point of view, one would expect that the SUSY is broken 
spontaneously so that the underlying microscopic Lagrangian is SUSY invariant but the vacuum state is not. This 
is the same mechanism which keeps the electroweak symmetry hidden at low energies. There is no consensus about 
how SUSY must be broken at high energies but by whichever way it is broken, at the electroweak scale the effect of 
SUSY breaking can always be parameterized by introducing additional terms in the Lagrangian which break SUSY 
explicitly. But these terms should only include operators which do not bring back the quadratic divergences from the 
quantum corrections to the scalar masses. These are called the SUSY breaking soft terms, all of which have positive 
mass dimension. 

Unfortunately, spontaneous breaking of SUSY in a phenomenologically acceptable way has not been achieved till 
now with only tree level renormalizable interactions. This problem has been solved by breaking SUSY in a sector, the 
so called "hidden sector", which has either no coupling or very small direct coupling to the MSSM fields. However, 
there are some mediating interactions which carry the information of SUSY breaking from the hidden sector to the 
MSSM and generate the soft terms. The gravitational interaction as a mediator of SUSY breaking has been the most 
popular scenario historically [6]. Another very popular mechanism has been the so called "Gauge Mediated SUSY 
Breaking (GMSB) " [7-12] scenario where the soft terms are generated through loops involving some messenger fields. 
These messenger fields couples to the SUSY breaking sector and also have SM gauge group quantum numbers. As 
gauge interactions are flavor diagonal, this mechanism also has the advantage of automatically being flavor-blind 
which is phenomenologically attractive from the points of view of Flavor Changing Neutral Currents (FCNC). 

The General gauge mediation (GGM), a framework developed by Meade, Seiberg, and Shih [13] generalizes models of 
gauge mediation to accommodate arbitrary hidden sectors including those which are not necessarily weakly coupled. 
It was shown in GGM that the soft masses of MSSM gauginos and sfermions, to leading order in the SM gauge 
interactions, can be expressed in terms of the hidden sector current correlation functions. Extension of this framework 
has been studied in [14-20]. It was also shown in a model independent way that non-universal gaugino masses can 
be obtained in GGM without spoiling the gauge-coupling unification. By universality in gaugino masses we mean 
that these masses at any scale are proportional to the corresponding gauge couplings at that scale. However, in 
weakly coupled GMSB scenarios with complete multiplets of SU(5) messengers the gaugino masses remain universal. 
In general messenger models [21, 22], one can achieve non- universal gaugino masses but the masses of the messenger 
fields get constrained in order to maintain gauge coupling unification. 

It was shown in [23] that one can achieve unification of gauge couplings even if the new fields added at some 
intermediate scales do not form complete multiplets of SU(5) and the unification is independent of the scale at which 
these fields are added. These specially chosen fields were termed as "Magic" fields and they were used in ordinary 
gauge mediation in order to get non universal gaugino masses with extremely large hierarchy, 1:30:200 among the 
U(l), SU(2) and SU(3) gaugino mass parameters [23], which is phenomenologically not interesting. Also the obtained 
patterns of mass ratios were not very flexible. 

In this paper we generalize their idea and write the supcrpotcntial containing all possible renormalizable terms for 
the messenger sector. We derive the soft SUSY breaking terms and show that the gaugino sector can be described 
by two independent parameters which can be tuned to achieve almost any ratio of gaugino masses. This is a novel 
feature of our models. We also show that the sfermion sector in this class of models can also be parameterized by two 
parameters which can be chosen to adjust the relative masses of squarks and sfermions. In this class of models the 
squark masses can be equal to the Higgs mass parameters rn Hu and m H at the messenger scale which can lead to a 
small value of the Higgsino mass parameter fi and thus, a Higgsino-likc NLSP (like models of extra ordinary gauge 
mediation (EOGM) with doublet triplet splitting [24]). Some region of the parameter space of this model also allows 
for the lightest SUSY Higgs mass about 125 GeV. 

We consider a phenomenologically interesting scenario where the ratio of the U(l) and SU(2) gaugino mass pa- 
rameters at the electroweak scale is comparatively larger than the universal case 1:2. In this case, because of the 
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comparatively larger splitting between the next-to- lightest neutralino lightest an d the neutralino (x?), the Higgs 
boson (h) coming from the decay \\ ~ > Xi h is quite boosted. Motivated by this, we study the £ + bb + fx channel 
to look for SUSY Higgs boson using the jet substructure technique. 

The paper is organized as follows. In the next section we briefly review the concept of magic fields. In Sec. Ill, we 
describe our models and then study the mass spectrum of gauginos and sfermions. The phenomenology of an explicit 
model of this class of models is discussed in Sec. IV A. We close our discussion in Sec. V with a brief summary of 
our findings. 



II. MAGIC FIELDS AND UNIFICATION 



In this section we will review the concept of magic fields and derive some relations which will be useful later in this 
paper. It is well known that, if no additional matter field is added in the intermediate scale, the SM gauge couplings 
get unified in the MSSM at the scale Mq ~ 2.0 x 10 16 GeV. In general, assuming that the gauge couplings unify at 
some scale Mqut, their one loop Renormalization Group (RG) equations reads, 

a- 1 (M) = -^ln— t—+ a - 1 (M GUT ), (1) 
Zn Mgut 

where the index a = 1,2,3 represent the SM gauge groups t/(l), SU(2) and SU(3) respectively and b a are the 
corresponding beta functions. In the MSSM we have b a = bg = {4p, 1, —3}. 

At any value of /x, Eq. (1) is an equation of a straight line in a -1 and b. It directly follows that for the unification 
of the gauge couplings, a^ 1 and b a must respect the following relation, 

c^ 1 - a^ 1 b 3 - b 2 ^ 



a — (Xi 



b 2 - b 1 ' 



Even if matter fields are added at some intermediate scale, the gauge coupling unification can be maintained (at one 
loop) so far as Eq. (2) is satisfied. The fields whose beta functions are such that they do not spoil Eq. (2) arc called 
Magic fields. 

Using the one-loop RG equation of Eq. (1) in the MSSM one has, 

a-\M z )-a^{M z ) = -M_j£ log ^£. (3 ) 

We now add matter fields at some intermediate scales and assume that fields added at the scale contribute to the 
beta functions by an amount df (dynkin index of that field). The gauge couplings at any scale \i can now be written 

as, 

a-\») = a-\M z )-f- log ±— £ flog A. (4) 

i with mKfi 

If we now choose 

d? - 4 = h (b a - b b ) , (5) 

then it is clear that Eq. (2) will remain intact. Hence we can achieve gauge coupling unification even with the- presence 
of intermediate mass scales and this is independent of the values of to,;. But unlike complete multiplets of say, SU(5) 
the unification scale changes in this case [25] . Let us assume that the gauge coupling constants now unify at a scale 
Mgut- We can now rewrite Eq. (3) as, 

a a (M z )-a b (Mz) = 7 — log— ^ 9 ^ lo § T7 ■ ( 6 ) 

27r Mqut . — zir Mgut 

i— I,-- - ,N 

Comparing Eq. (3) and (6) we can now write the new unification scale Mqut in terms of the unification scale Mo in 
the MSSM, 

Mgut = M J] ( ^ ) , where m = h/ I 1 + ^ kA . (7) 



t=l,2,.» .N 
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The change in the unified value of the gauge couplings, da 1 (Mqut) can also also be written down explicitly, 

7T i=l,2,— ,N\ j=l,2,— ,N J 

An example of magic fields is the combination of fields <f>Q , 4>q and 4>g whose transformation properties under the 
SM gauge group are given by (3,2)i, (3, 2)_i and (8, l)o respectively. Here the first number in the bracket is the 

6 6 

SU(3) representation and the second number refers to the SU(2) representation. The number outside the bracket is 
the Hypercharge of the multiplct. Note that the individual fields above are not magic fields but the combination as 
a whole satisfies the magic field condition of Eq. (5). The set of fields mentioned above can be obtained from the 
SO(10) complete multiplcts 16, 16 and 45 after spontaneous symmetry breaking [23]. 

III. MASS SPECTRUM OF THE MAGIC MESSENGER MODELS 

In this section, we will derive some general properties of the magic messenger models where, unlike models of 
ordinary gauge mediation, the supcrpotcntial also contains bare mass terms allowed by symmetry. But before we go 
to the details of our model, let us first discuss those models of GMSB where non-universal gaugino masses have been 
obtained and their qualitative differences with our model. 

In ordinary gauge mediation models the messenger fields are complete multiples of SU(5) like 5 © 5. Addition of 
such complete multiplets changes the beta functions equally for all the SM gauge groups at any scale. Clearly this 
does not change the universality of the gaugino masses. Even if all possible renormalizable terms are considered in 
the superpotential, the situation does not change and universality is maintained. But instead of SU(5) multiplcts 
if one uses SU(2) and SU(3) irreducible representations and write the most general superpotential (like models of 
EOGM with doublet triplet splitting [24]), then non- universal gaugino masses can be achieved only if these doublets 
and triplets arc charged differently under some global symmetry. But to achieve gauge coupling unification at the 
same time the masses of the messenger fields should satisfy some constraints. In generalized messenger models, one 
can easily obtain non-universal gaugino masses but this type of models generally bypass the question of unification 
by arguing that there may be some other fields (not messenger of SUSY breaking) in the theory which are charged 
under the SM and using these fields it is always possible to get successful unification. If one does not assume this 
then messenger masses get constrained again. Another generalized messenger model is also available in the literature 
[26] where messenger masses are not constrained but the GUT used there is an anomalous U(l) GUT. 

The models we are proposing are also generalized messenger models but here the messenger fields are magic fields. 
This class of models has two main differences with other generalized messenger models: (a) gauge coupling unification 
is independent of messenger scales and (b) changes in beta functions for each added messenger satisfy the magic 
relation of Eq. (5). The second property has non-trivial consequences as we will see in the next two sections. 

A. The models 

Our models consists of (1) ./V pairs of magic messenger fields, each pair consisting of ((/>, <f>) where cj) {<t>) transforms 
under some representation (its conjugate representation) of the SM gauge group, and (2) a spurion field X which is 
a standard model gauge singlet but can have charges under some other symmetries like the R-symmetry. Note that 
these 4> fields are in general not single irreducible representations; they can be a set of fields (each transforming under 
an irreducible representation) which satify the magic condition of Eq. (5). We assume that the field X gets a vacuum 
expectation value (VEV) through some dynamics of the hidden sector: (X) = M + 9 2 F. For the sake of calculational 
simplicity we make the assumption that two different pairs of magic messengers either transform under the same 
representation of the gauge group or have no constituent fields having the same transformation property. This allows 
us to sub-divide the messenger field pairs into different sets (labeled by p) where all members of a set have the same 
representation. Note that a set may consist of either one pair of messenger fields or many pairs of them. We label the 
different members (a particular pair of magic messengers) of a set p by the index i p . For example, if a particular set 
p has five pairs of magic messengers then i p runs from 1 to 5. We now write the renormalizable superpotential to be, 

P ip,j P 

where the form of the matrices X p and m v are determined by global (R and/or non-R) symmetries of the theory. 
Note that, in general different constituent fields of any magic field can get different masses and couplings. But to 
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form a magic field, these constituent fields need to share the same mass. We treat the magic fields (fa 's) 'as if they 
are irreducible representations of some group' and write the couplings in the superpotential because this assumption 
automatically ensures the above requirement (of equality of the masses of the constituent fields of each the magic 
messenger fields). The above superpotential is symmetric under the interchange of fa p field with fa p and vice versa. 
This symmetry is called the messenger parity and it helps us to get rid of the dangerous Fayet-Iliopoulos (FI) terms. 
1 Throughout this paper we will call this class of models as "Magic Messengers in Gauge Mediation (MMGM)". 

To calculate the soft-terms, we need to write the Lagrangian in mass eigenstates of fermion fields and scalar fields. 
Using bi-unitary transformations on the superfields, we first make fermion mass matrices of each set diagonal and 
real: Mj? = diag(- • • , m ip , ■ ■ ■). In this basis the Kahler term will not change. The matrices X P F will change but for 
brevity we keep the same symbols. Messenger parity and CP conservation imply that these matrices should be real 
and symmetric. Going to the basis </>± p — ^75 (' ' • 1 4>i P i 0* p j ' ' ' ) T j we can bring sfermion mass squared matrix of any 
set in the block-diagonal form with two blocks: M± p = (Mj? ) 2 ± X P F. Now these two matrices can be diagonalized 

as U± M± p U± with eigenvalues . We define the following combinations of U matrices which will be used in the 
next section, 

= u ±U uP ±^ and - E(KikU- kpjp )(u%u +lpip ) . (10) 



B. Gaugino masses 



The expression for gaugino masses at the messenger scale can be written as, 



p 



(11) 



where is the dynkin index of the magic messenger pairs in the set labeled by p corresponding to the gauge group 
label a, the superscript G in A refers to the gaugino sector, and A p is given by [22] 



Ag= 2 v(±Mii ^ m± i- iog 



m 2 ±Jp 



(12) 



Here A p are independent of the SM gauge group. However the presence of dp in the expression for gaugino masses 
implies that the ratio of gaugino masses is not equal to the ratio of a Q 's. This means that the gaugino masses are 
non-universal in this class of models. Apparently it seems that the gaugino masses are completely independent, but 
actually this is not true. Their ratio satisfies some beautiful structure. The ratio of the gaugino masses can be written 
as, 



Mi : M 2 : M 3 = 1 



a 2 



1 



(2) 



where 



c 



«3 



1 



,0) 



The value of C can be zero, positive or negative. All the gaugino mass parameters become equal if 



(13) 



(14) 



C = Co 



a 3{0 -t>o ) ~ a 2{0 -O ) 



«2(4 2) 



where in the last step we used Eq. (2). Any arbitrary £ can always be written as 

C = Co + c 



(15) 



1 This FI term is proportional to the VEV of the scalar component of the U(l) current superficld J. Under messenger parity J — > — J, 
because the fields <f>i p transform under the conjugate representation of the fields <j>i p . Hence, due to messenger parity, VEV of J vanishes 
[13, 16]. 
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In that case Eq. (13) takes the form, 



M 1 :M 2 :M 3 = l:l + ^^-4 1) )C:l + ^^ 3) -^ ) )C, 



a 5 o;i 



It is now clear from Eq. (16) that various ratios among the gaugino masses can be obtained, (a) for negative Q, one 

48q 3 : 



gets normal hierarchy, (b) for £ = 0, all the gaugino masses are equal, and, (c) for non-zero positive C < one 



gets inverted hierarchy. For £ > some more patterns can be achieved. 

So to conclude this section, the gaugino sector of this class of models can be parametrized by two free parameters: 
one of them can be taken as the mass of the U(l) gaugino and the other one is the £ (or £) parameter. These two 
parameters can be tuned to get phcnomcnologically interesting ratios among the gaugino masses. 



C. Sfermion masses 



The sfermion masses are obtained at the two loop level in GMSB. In GGM, these are parameterized by three 
parameters A a (a =1, 2, 3) for three gauge groups U(l), SU(2) and SU(3). In terms of these three parameters the 
expression for the sfermion masses [13] can be written as follows, 
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g c 2 (f,a)A a , 



(17) 



where C2(/, a) is the quadratic Casimir of the representation / of the gauge group a. 
expression for A a turns out to be 



In our class of models the 



A„ 



1 



(16tt 2 ) : 



E d P A P 



(18) 



where is given by [22], 



A 



2 E m U 

ip jp =t 



—Bi i Li 2 1 

2 % VJV z i 



Jp 



-2A P A Li 2 1 



tpjp 



(19) 



The quantities A^ and A a are required to be strictly positive in order to have positive sfermion masses. 
Ratio of the A a parameters can be written as, 



28 48 
At : A 2 : A 3 = 1 : 1 - — r\ : 1 - — r], 
5 5 



(20) 



where 



(21) 



Note that 77 cannot be greater than rjo = j^. For 77 ~ 770, one has A\ > A 2 ^> A3. In this limit, from Eq. (17) one can 
see that soft SUSY breaking squark and Higgs masses will be approximately same (universal) at the messenger scale. 
This would in turn render the value of m H at weak scale insensitive to the universal scalar masses in the ultraviolet. 
This can be understood as a consequence of a "focus point" in the RG behavior of m H [27, 28]. In this case it is 
possible to have a small fi (supersymmetric mass term for the Higgs) [29, 30] and thus, a Higgsino like NLSP. This 
can lead to very distinct phenomenology [31, 32]. 
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IV. AN EXPLICIT MODEL OF MMGM 



In this section we construct an explicit model of MMGM. We take two sets of messengers, the first set constists of 
two pairs of SU(5) 5 + 5 fields: <pi, 4>i with i = 1,2, and the second set contains a magic field: 3 = tpq + <j>Q + <pQ. 
Note that one can write down the required terms for gauge mediation (as given in Eq. (9)) using only </>3 and without 
invoking its conjugate field fa. Moreover, absence of the field </>3 does not break the messenger parity because the 
superpotial remains invariant under the interchange of 4>q by </>q and vice versa. This is why we do not include the 

(j>3 field in the second set just to make it more economical. In the second set, we will denote the fermion mass by 
m 3 and the eigenvalues of the scalar mass squared matrix by m§ ± d. For the first set, in the diagonal basis of the 
fermion mass matrix Mf = diag(mi, TO2), the mass squared matrix of the scalars will look like, 

Clearly the eigenvalues and the diagonalizing matrices of M± will involve the quantities a, b, c and hence, the four 

parameters Mi, £, A\ and rj will also depend on them. Just to give few examples, in Table-I we choose six benchmark 
points 2 and show the values of a, b, c, d and the corresponding numerical values for Mi, A\ and 77. 



Benchmark 
points 


a 


b 


c 


d 


Mi 
(GeV) 


C 


(10 s GeV 2 ) 


V 

(10- 8 ) 


1 


0.85 


80 


0.85 


0.60 


97.304 


-0.1914 


2.0534 


-1.4061 


2 


0.80 


70 


0.80 


0.90 


95.165 


-0.2794 


1.5722 


-4.1320 


3 


0.80 


70 


0.80 


1.2 


98.373 


-0.3527 


1.5722 


-7.3458 


4 


0.80 


70 


0.80 


1.4 


100.51 


-0.3989 


1.5722 


-9.9983 


5 


0.80 


70 


0.80 


1.6 


102.65 


-0.4432 


1.5723 


-13.059 


6 


0.80 


70 


0.80 


1.8 


104.79 


-0.4858 


1.5723 


-16.527 



TABLE I. The values of a, b, c and d parameters for six benchmark points and the corresponding numerical values for the 
quantities Mi, C, Ai and rj. The values of mi, m,2 and 7713 are 1.0 x 10 14 GeV. The parameters a, b, c and d are given in units 
of 10 18 GeV 2 . 

It is well known that among the various models of GMSB, the mass ratio between the sfermions and gauginos is 
the lowest in case of ordinary gauge mediation models [20] (which can be obtained by setting b = in Eq. (22)). In 
order to increase this ratio, b should be made large compared to a and c. The parameter d is related to the magic 
part of the model. So changing d, one can increase or decrease the splitting among the gaugino masses. Note that, we 
got ( ~ 1CP 1 whereas 77 ~ 1CP 5 . So gaugino masses are highly hierarchical (see Table II) but A a parameters are not. 
To understand this one should note that the SU(5) part of this model is not an OGM model but an EOGM model, 
and b is hundred times larger than a, c and d. From Eqn. (B.6) and (B.8) of [20], one can see that the off-diagonal 
element b has no effect on the expression for gaugino masses whereas it has dominant contribution on the sfermion 
masses. Using the fact that ~ 10~ 8 , the expression for the 77 parameter can be approximated as: 77 ~ — This 
explains why 77 is small in this model. 

A. Phenomenology: boosted higgs signal 

There are many models of supersymmetry breaking where the soft supersymmetry breaking gaugino mass terms 
Mi, M2 and M3 meet to a common value m.1/2 at the GUT scale. Now the one-loop Renormalization Group equations 
for the three gaugino mass parameters in the MSSM are determined by the same quantities b^ (a=l, 2, 3) which also 
control the RG running of the three gauge couplings. It then immediately follows that each of the three ratios ^ 
is one-loop RG invariant. In models of gaugino mass unification (for example, models with minimal supergravity or 
gauge-mediated boundary conditions) this leads to an interesting relation, Mi:M 2 :M 3 = 1:2:7 approximately at the 
TeV scale (modulo two loop corrections and unknown threshold corrections). Now, if the supersymmetric Higgsino 



2 In our models, gravitino mass is in the range 17 GeV to 20 GeV. If we assume gravitino to be the LSP then Big Bang Neucleosynthesis 
will be problematic because of large neutralino lifetime. This problem can be easily solved in the axino dark matter models where life 
time of the lightest neutralino can be less than .1 sec [33-35]. 
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mass parameter /i 3> Mi, M2 then the physical mass eigenstates consist of a "bino-likc" lightest neutralino Xi and a 
"wino-like" next-to-lighest neutralino x° an( i lightest chargino xt . In this case one has to a very good approximation 
nig ~ M3, m^o,m~± ~ M2 and m^a ~ Mi. With the increase of the lower bound on squark and gluino masses by 
the LHC data, their production cross section has been pushed to quite low values (~ few fb). On the other hand the 
electro weak gaiginos can be sufficiently lighter (xt > X2 ~ 250 is still allowed even if universality is assumed) . Hence 
the pair production cross-section of the light electroweak gauginos dominate the SUSY production cross-section. 

In our model, as we have already seen, in general the gaugino mass parameters can be arranged to have any ratio 
among themselves. In particular, we consider the case where Mi, M2 and M3 are much more hierarchical than the 
ratio 1:2:7. In this case Xi an d X2 even lighter than 250 GeV are still allowed. In this section we consider the 
production of X11X2 ~ an d their subsequent decays, xt ~^ XiW and X2 ~~ ^ Xi^ 1 - The prospect of this channel at 
8 TcV LHC in case of the mSUGRA model has been studied in detail in [37]. It was concluded that an integrated 
luminosity of 100 fb -1 will be needed for a good signal to background ratio. Note that, unlike mSUGRA in our model 
the lightest neutralino Xi can be considerably lighter than x\ an d hence, the lightest Higgs boson from the decay of 
X° can be quite boosted. 

In Table-II we show a few benchmark points which will be used for a detail signal anlysis. We choose the parameters 
in our model so as to get different ratios of gaugino masses along with the lightest Higgs boson mass consistent with 
the recent hints of Higgs signal by the CMS [39] and ATLAS [40] collaborations. Note that, in our case the branching 
ratio for the decay X2 ~~ ► Xih is very large (> 95%). For related discussions on the parameter dependence on this 
branching ratio and interplay with other decay modes see the reference [41]. 

To generate the mass spectrum of the SUSY particles we have used the package SuSpect [42] . SUSYHIT [43] has 
been used to calculate the corresponding branching ratios. Note that the gluinos are quite heavy for all the benchmark 
points and consequently their production cross-section is extremely small. For example, for the benchmark point-1 
the gluino pair production cros-section is about 4 fb (LO) for 8 TeV center of mass energy and it is even smaller for 
the other benchmark points. 



Benchmark 
points 


tan/3 


m~± 

Xl 

(GeV) 


m-o 
(GeV) 


m-o 
(GeV) 


m h 
(GeV) 


m- g 
(GeV) 


(GeV) 


BR(x^ -»■ X?ft) 




(10 4 ) 


1 


20 


203.0 


203.0 


50.9 


124.8 


1019 


7935 


98.908 


66.25 


1.53 


2 


25 


248.1 


248.1 


49.0 


124.6 


1239 


7003 


96.883 


30.83 


1.19 


3 


30 


299.4 


299.4 


50.2 


124.8 


1482 


7024 


95.221 


19.74 


1.19 


4 


20 


332.3 


332.3 


50.6 


124.8 


1641 


7055 


96.076 


52.64 


1.21 


5 


20 


366.1 


366.1 


51.3 


124.9 


1797 


7070 


95.621 


52.92 


1.21 


6 


20 


399.7 


399.7 


51.9 


124.9 


1950 


7085 


96.343 


53.19 


1.22 



TABLE II. Masses of the lightest neutralino, next-to-lightest neutralino, lightest chargino and lightest CP even neutral Higgs 
for the six benchmark points. All other SUSY particles have masses in the multi-TeV range (~ 7 TeV to 15 TeV). Definition 

of fine tuning parameters are ^M = §[l + fr ff'^I'^ ] ^ and ^ = 4 ^ ^ P m^IT}-^ where *P = 
(tan 2 P + l)/(tan 2 ft — 1). Hence fine-tuning is large. 

In Fig.l we show the transverse momentum distribution of the Higgs in the decay X2 ~~ ^ Xi^ 1 following the direct 
production of Xi , X2 & t the LHC with 8 TeV center of mass energy. It can be observed that a large fraction of the 
Higgs bosons has transverse momentum greater than 100 GeV. This allows us to use the jet substructure technique 
to look for Higgs in the decays of directly produced electroweak gauginos. 

The use of jet substructure for the reconstruction of hadronic decays of boosted W, Z, Higgs bosons and top quarks 
has received considerable attention in recent years [45-47]. A study of jet substructure in the context of a search 
for a heavy Higgs boson decaying to W W was first carried out in Ref. [48]. More recently, Butterworth, Davison, 
Rubin and Salam(BDRS)[49] studied the case of a light Higgs boson (rriH ~ 120 GeV) produced in association with 
an electroweak gauge boson. The leptonic decay of the associated vector boson provides an efficient trigger for these 
events. The BDRS algorithm involves a technique using the mass-drop and the filtering to transform the high-p T 
WH, ZH(H -> b b) channel into one of the best channels for discovery of Standard Model Higgs with small mass at 
the LHC. 

In this section we adopt the BDRS method for tagging hadronically decaying Higgs boson. We describe below the 
exact procedure adopted in our analysis in order to implement this along with our other selection cuts. 

We first cluster all the stable final state particles (excluding leptons, neutrinos and neutralinos) into "fat-jets" 
using the Cambridge-Aachen algorithm (CA algorithm) [50, 51] as implemented by the Fastjet package [52] with R 
parameter of 1.2. We then select the jets with transverse momentum p T > 100 GeV and pseudo rapidity \rj\ < 2.5. 
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FIG. 1. Left panel : The transverse momentum distribution of the Higgs boson from the decay \2 ~ ► Xi h following the direct 
production of (xt,X°) at 8 TeV LHC for our SUSY benchmark point 5. The y-axis has been normalized to 1. 
Right panel: The invariant mass distribution of the reconstructed Higgs boson using the jet substructure algorithm, as described 
in the text, for our SUSY benchmark point 5. The y-axis has been normalized to 1. The figures have been generated using the 
CERN package PAW [44]. 



We now perform the jet substructure analysis on the hardest jet following the BDRS prescription. Here we closely 
follow the discussion of Rcf. [49], mentioning our specific choices of parameters as and when the occasion arises. The 
first step is to select a jet-j and apply the following procedure: 

1. Undo the last clustering step of the jet-j to get two subjets. Label the two subjets ji and j2 so that > 
mj 2 . (Remember that at each step during the CA clustering the masses of the proto-jets to be combined are 
recorded.) 

m.j min(p^ ■ , p^, ■ ) 

2. Define /x= — — , y cu t= ^ ^^ji h' Here ARjj i j 2 is the separation between the two subjets ji 

mj im 

and j2 in the pseudo rapidity^)- azimuthal angle(0) plane. If fi < fx cu t (significant mass drop) and y > y cu t 
(the splitting of the hard jet-j into two subjets ji and j2 is not too asymmetric) then go to step 4. We use ^ cut 
=0.4 and y cu t=0.1 in our analysis. 

3. Otherwise redefine j = ji and go back to step 1. 

4. Take the constituents of the mother jet-j and recluster them with the Cambridge-Aachen algorithm with an 

R-parameter of Rfiit=miri( Aflj 2 1 ' 12 ,0.4). Construct n new subjets jf", jf 1 *, jf 1 *, jJJ 1 * ordered in descending 

Pt- This step is supposed to reduce the degradation of resolution on jets caused by underlying events. 

5. Require the two hardest of the subjets jf 1 * j^ 1 ' to have b tags. 

min(n,3) 

6. Define j hl ss s = j? 1 *. This step captures the dominant 0(a s ) radiation from the Higgs decay, while 

i=l 

eliminating much of the contamination from underlying events [49]. 

The parameters involved in this method can be, in principle, optimized event by event [53]. In our analysis, 
these parameters have been set to fixed values as already mentioned above. In Fig. 1 we show the distribution of 
the reconstructed Higgs mass using the above BDRS prescription for our SUSY benchmark point-5. For the other 
benchmark points the distributions are qualitatively same and we do not show them here. 

We also impose a b-jet reconstruction efficiency of 70% [54] in our analysis. We then require the following pre- 
selection cuts [36-38]: 
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Number of Events after individual cuts 




Process 


Crosssection 
(LO) 


Simulated 
Events 


Cut-I 


Cut-II 


Cut-III 


Cut-IV 


Cut-V 


Crosssection 
after Cut-V 


K 

—=(30 fb ) 


Signal 




Benchmark- 1 


533 fb 


50K 


681 


178 


71 


47 


28 


0.21 fb 


ss 4 


Benchmark-2 


241 fb 


50K 


1220 


311 


150 


106 


84 


0.27 fb 


« 5.1 


Benchmark-3 


108 fb 


50K 


1729 


475 


287 


201 


177 


0.27 fb 


« 5.1 


Benchmark-4 


71 fb 


50K 


2112 


552 


356 


270 


246 


0.23 fb 


« 4.4 


Benchmark-5 


45 fb 


50K 


2697 


788 


522 


380 


341 


0.21 fb 


« 4 


Benchmark-6 


30 fb 


50K 


3092 


904 


655 


488 


449 


0.21 fb 


« 4 


Background 




t t(0-100) 


48.3 pb 


40L 


1126 


274 


29 


3 







t t( 100-200) 


36.6 pb 


30L 


976 


234 


26 


6 


1 


0.0119 fb 


t t (200-300) 


7.8 pb 


5L 


157 


40 


6 


3 


2 


0.0306 fb 


t t (300-500) 


1.8 pb 


1.5L 


24 


5 


1 










t t(500-oo) 


134 fb 


10K 



















W(-> lv,) b b 
I = e, /Li, r 


3 pb 


201308 


64 


35 


3 


2 







W h 


549 fb 


50K 


401 


125 


9 


5 


2 


0.0108 


W Z 


13 pb 


8L 


11 


3 













Z h 


296 fb 


20K 


167 


4 













z(-> rr)bb 

l = e,fi,r 


2 pb 


126581 


26 


5 













t(— > e ^ e b) b 


308 fb 


36817 


6 


5 













t b W 


18.7 pb 


597812 


156 


43 


4 


2 


2 


0.0306 fb 


Total 
Background 




0.084 fb 



TABLE III. Event summary for the signal and the backgrounds after individual cuts as described in the text for LHC8. In the 
second column, the Leading Order (LO) cross-sections have been obtained using PYTHIA and ALPGEN. While calculating 
the final crosssection after Cut-V, the NLO cross-sections for signal(from PROSPINO) and appropriate K factors (whenever 
available) for the backgrounds, as mentioned in the text, have been used. The b-tagging efficiency has also been multiplied. K 
and L in the third column stand for 10 3 and 10 s respectively. Note that the number of simulated events is always more than 
the expected number of events at LHC (at 30 fb -1 ) for both the signal and backgrounds. 



• Cut-I : We require the mass of j hl ss s to be in the window [119, 129]. 

• Cut-II : Exactly one isolated lepton(£) with pt(^) > 20GeV and no isolated lepton with 10 GeV < pt(^) < 20 
GeV. 

• Cut-III : The transverse mass of the lepton(^)-missing Pt system M^ T ' > 90 GeV. 

• Cut-IV : At this stage we again construct "normal" jets using the CA algorithm with R parameter of 0.5, |?y| < 
2.5, Pt > 50 GeV. We then calculate Ht, which is defined as the scalar sum of the Pt's of all these "normal" 

jets. We define R T b = ^ Tbl P Tb2 . Remember that pTt>i and PTb 2 are the transverse momenta of the two 

subjets j^ 1 * and j| which are by now identified as two b-jets. We demand R T b > 0.9. 

• Cut-V : Events are selected with Pt> 125 GeV. 

In Table-Ill and Table-IV we show the signal and backgrounds after each selection cut for 8TeV and 14TeV center 
of mass energies respectively. For all the signal points we use the NLO cross-section from Prospino[55]. We simulate 
the tt, Wh, WZ and Zh backgrounds using PYTHIA[56]. For the Wbb, Zbb and single top backgrounds we generate 
the unweighted event files in ALPGEN [57] and then use the ALPGEN-PYTHIA interface (including matching of 
the matrix element hard partons and shower generated jets, following the MLM prescription[58]) to perform the 
showering and implement our event selection cuts. For the tt background a K factor of 2 has been used. We use the 
CTEQ6L parton distribution functions and set top mass at 172.9 GeV in our analysis. We have also used the new 
LHC PYTHIA 6.4 Tune Z2* for the correct description of the Underlying events[59]. We see that a S/VB about 
4-6 can be obtained at 8TeV LHC with an integrated luminosity of about 30 fb -1 . At the 14TeV LHC the siuation 
is much better and even with 15 fb -1 we get fairly good number of events, while the backgrounds are totally under 
control with our selection cuts. In Table-Ill (Table-IV) we have simulated at least 30 fb" 1 (15 fb -1 ) of events for 
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Number of Events after individual cuts 




Process 


Crosssection 
(LO) 


Simulated 
Events 


Cut-I 


Cut-II 


Cut-III 


Cut-IV 


Cut-V 


Crosssection 
after Cut-V 


Signal 


Benchmark- 1 


1.4 pb 


50K 


769 


186 


78 


54 


29 


0.53 fb 


Benchmark-2 


681 fb 


50K 


1282 


327 


156 


119 


92 


0.73 fb 


Benchmark-3 


332 fb 


50K 


1857 


465 


288 


198 


169 


0.73 fb 


Benchmark-4 


226 fb 


50K 


2363 


626 


392 


274 


238 


0.73 fb 


Benchmark-5 


155 fb 


50K 


2792 


746 


499 


359 


321 


0.67 fb 


Benchmark-6 


105 fb 


50K 


3125 


843 


592 


416 


373 


0.53 fb 


Background 


t t(50-oo) 


335 pb 


100L 


3467 
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11 


3 


0.098 fb 


W(-> I v t ) b b 
I = e, fx,r 


5.45 pb 


429871 
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86 


7 


6 







W h 


1.24 pb 


50K 
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7 


3 


1 


0.012 fb 


W Z 


29 pb 


8L 


20 


5 












Z h 


674 fb 


50K 


449 


6 













Z(-» r^~)bb 
t = e,fi,T 


7.4 pb 


607607 


116 


23 


5 


4 







t b 


5.6 pb 


182974 


16 


3 













t(-> ^b) b 
W (-> had) 
£ = e,fj,,r 


17.1 pb 


499383 


125 


65 


11 










t(-> had) b W 
(->^) 
£ = e, /i, t 


17.1 pb 


665454 


229 


127 


14 


3 







Total 
Background 




0.11 fb 



TABLE IV. Event summary for the signal and the backgrounds after individual cuts as described in the text for LHC14. In the 
second column, the Leading Order (LO) cross-sections have been obtained using PYTHIA and ALPGEN. While calculating 
the final crosssection after Cut-V, the NLO cross-sections for signal(from PROSPINO) and appropriate K factors (whenever 
available) for the backgrounds, as mentioned in the text, have been used. The b-tagging efficiency has also been multiplied. K 
and L in the third column stand for 10 3 and 10 s respectively. Note that the number of simulated events is always more than 
the expected number of events at LHC (at 15 fb -1 ) for both the signal and backgrounds. 

both the signal points as well as the backgrounds. Note that detector and other experimental effects will degrade the 
signal significance somewhat. We have checked that the numbers change by ~10% if a Gaussian smearing is added 
to the transverse momenta of the jets and the missing transverse momentum, though a faithful quantification of the 
detector effects is beyond the scope of this work. 

As the masses of Xi and x° increase (see Table. II), their production cross-section gradually decreases. This tends 
to reduce the signal to background ratio. But on the other hand from the benchmark-1 to benchxnark-6 the mass 
difference between \2 an d Xi a l so increases. This makes the Higgs boson more boosted thereby increasing the efficiency 
of the jet substructure algorithm. With the increasing mass of x® i n Table. II these two opposite effects keep competing 
with each other. This is why the signal efficiency initially increases with increasing x_2 m &ss but again starts falling 
down because of the rapid decrease in the crosssection. 

V. CONCLUSION 

In this paper we have implemented magic fields as messengers of SUSY breaking in GMSB. One of the advantages 
of using magic fields as messengers over other generalized messengers is that achievement of unification is independent 
of the masses of the magic messengers. 

In our model the gaugino sector is parameterized by only two independent parameters, one of them can be taken 
to be the U(l) gaugino mass and the second one being the ( parameter (Eq. (14) and (15)). The £ parameter can be 
tuned to get various hierarchies among the gaugino masses which can lead to distinct phenomenological consequences. 

The sfcrmion sector can also be characterized by only two independent quantities. These arc the U(l) A-parameter 
A\ (Eq. 18) and the parameter r\ (Eq. 21). Choosing 77, different hierarchies between the squark and slepton masses 
can be achieved. When the value of rj is close to its upper limit 770, the squark and Higgs masses at the messenger scale 
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tend to be almost same. This allows one to have small ji parameter and Higgsino-like NLSP similar to the models of 
EOGM with large doublet-triple splitting. 

We focus on the region of parameter space where a comparatively larger splitting (about 1:6) between the U(l) 
and SU(2) gaugino masses is achieved along with the lightest supersymmetric Higgs boson mass about 125 GeV. We 
consider the direct electroweak production of xt and x\ with xt decaying to the lightest neutralino Xi and a W 
boson, and x\ decaying to x\ and the lightest Higgs h. Because of the large splitting between x° and Xii the produced 
Higgs boson is expected to have quite large transverse momentum. Motivated by this we have analyzed the I + b b + 
fx channel using the jet substructure technique. We have simulated all possible backgrounds for this final state and 
conclude that while S/V~B ~ 4 - 6 is viable (for the mass ranges of charginos and neutralinos we have considered) 
at 8 TcV LHC with a integrated luminosity of 30 fb _1 , LHC14 can do much better and even with 15 fb _1 of data a 
decent number of signal events over the backgrounds is expected. In our analysis, we have not considered any detector 
effect which is expected to degrade the signal significance to some extent. 

A detailed exploration of other phenomenological consequences of this class of models including constraints from 
Flavor physics as well as other low energy experiments should be carried out and we plan to perform such a study in 
a future publication [60]. 
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